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We study Hamiltonian systems with point interactions and give a systematic descrip- 
tion of the corresponding boundary conditions and the spectrum properties for self-adjoint, 
PT-symmetric systems and systems with real spectra. The integrability of one dimensional 
many body systems with these kinds of point (contact) interactions are investigated for 
both bosonic and fermionic statistics. 
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The complex generalization of conventional quantum mechanics has been in- 
vestigated extensively in recent years. In particular it is shown that the standard 
formulation of quantum mechanics in terms of Hermitian Hamiltonians is overly 
restrictive and a consistent physical theory of quantum mechanics can be built on 
a complex Hamiltonian that is not Hermitian but satisfies the less restrictive and 
more physical condition of space-time reflection symmetry (PT symmetry) It 
is proven that if PT symmetry is not spontaneously broken, the dynamics of a 
non-Hermitian Hamiltonian system is still governed by unitary time evolution. A 
number of models with PT-symmetric and continuous interaction potentials have 
been constructed and studied i2j. In this article we study Hamiltonian systems 
with singular interaction potentials at a point. We give a systematic and complete 
description of the boundary conditions and the spectra properties for self-adjoint, 
PT-symmetric systems and systems with real spectra. We then study the integra- 
bility of one dimensional many body systems with these kinds of point interactions. 

1 Self-adjoint point interactions 

Self-adjoint quantum mechanical models describing a particle moving in a local 
singular potential concentrated at one or a discrete number of points have been 
extensively discussed in the literature, see e.g. |5|^^5 and references therein. One 
dimensional problems with contact interactions at, say, the origin {x = 0) can be 
characterized by separated or nonseparated boundary conditions imposed on the 
wave function at x = The first model of this type with the pairwise inter- 
actions determined by (5- functions was suggested and investigated in |7]. Intensive 
studies of this model applied to statistical mechanics (particles having boson or 
fermion statistics) are given in [Sj |^ . 
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Nonseparated boundary conditions correspond to the cases where the perturbed 
operator is equal to the orthogonal sum of two self-adjoint operators in L2(— oo,0] 
and L2[0, cxd). The family of point interactions for the one dimensional Schrodinger 
operator can be described by unitary 2x2 matrices via von Neumann formulas 
for self-adjoint extensions of symmetric operators. The non-separated boundary 
conditions describing the self-adjoint extensions have the following form 

/ V(+0) \_,^ef a b\f V(-0) A ... 

[ v'(+o) j-^ [c d )[ v'(-o) ; ' 

where ad — be = 1, 9^a,b,e,d d H, 'ip{x) is the wave function of a spinless particle 
with coordinate x. The values 9 — b — 0,a~d^linQ correspond to the case of 
a positive (resp. negative) 5-function potential for c > (resp. c < 0). For general 
a, 6, c and d, the properties of the corresponding Hamiltonian systems have been 
studied in detail, see e.g. |Bl [TI11[T^ . 

The separated self-adjoint boundary conditions are described by 

^'(+0) = h+^P{+0) , ^'(-0) = h-iPi-O), (2) 

where E IRU {oo}. h+ = CX3 or /i — oo correspond to Dirichlet boundary 
conditions and ^ or ^ correspond to Neumann boundary conditions. 



2 PT-symmetric point interactions 

An operator is said to be PT-symmetric if it commutes with the product opera- 
tor of the parity operator P and the time reversal operator T. It can be shown that 
the family of PT-symmetric second derivative operators with point interactions at 
the origin coincides with the set of restrictions of the second derivative operator to 
the domain of functions satisfying the boundary conditions at the origin 

*'"°'^^«r,t'-°n (3) 



for non-separated type, where 



^ , VT+bce'^ 



c vT+bc e-*"= 

the real parameters b > 0,c > —1/5^), 9,(p G [0, 27r); or corresponding to the 
separated type 

hoi^'i+O) =- /iie'V(+0), ho^l>'i+0) = -/iie-'V(-O) (4) 

with the real phase parameter 9 e [0, 27r) and the parameter h = (ho, hi) taken 
from the (real) projective space P^. 

^) If the parameter b is equal to zero, then the second inequahty should be neglected. 
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The spectrum of any PT-symmetric second derivative operator with point inter- 
actions at the origin consists of the branch [0, oo) of the absolutely continuous spec- 
trum and at most two (counting multiplicity) eigenvalues, which are real negative 
or are (complex) conjugated to each other. The eigenvalues corresponding to PT- 
symmetric eigenfunctions are real and negative. Every eigenfunction corresponding 
to any real eigenvalue can be chosen either PT-symmetric or -antisymmetric. 

The spectrum of the PT-symmetric second derivative operator with non-separated 
type point interaction at the origin is pure real if and only if the parameters ap- 
pearing in (PJ satisfy in addition at least one of the following conditions: 

be sin^ < cos^ 0; (5) 



6c sin ip > cos </> and cos<p > 0. (6) 

3 Point interactions with real spectra 

We consider further non-separated type boundary conditions at the origin lead- 
ing to second derivative operators with real spectrum. A general form of the bound- 
ary condition can be written as 



^^(+0) \_f a (3\f V(-0) 
^'(-fO) J [7 6 )[ V'(-O) 



(7) 



where a, /3, 7 and J G C We suppose that the matrix B — [ ^ ] appearing 



7 S 

in the boundary conditions is non degenerate (from GL{2, <C)). Again it is easy 
to prove that the operator has branch of absolutely continuous spectrum [0, cxd). To 
study its discrete spectrum we use the following Ansatz for the eigenfunction 

^(x) = { : ; ' > 0, 

corresponding to the energy X = k^. Substituting this function into the boundary 
conditions (0 we get the dispersion equation + ik{a -f (5) — 7 = 0. 

The set of coefficients a, /3, 7, and S satisfying the condition 3A;i,2 < can be 
parameterized by 8 real parameters and leads to operators with pure absolutely con- 
tinuous spectrum [0,oo). Pure imaginary solutions to the dispersion equation leads 
to nontrivial discrete spectrum. Set t — a + 6. The solutions are pure imaginary if 
and only if the following conditions are satisfied: 

r = te'^ /? = 6e'^ 7 = ce^', 4^< ^, (8) 

for /3 7^ 0, where t, b, c are real numbers. If /3 = 0, the spectrum is guaranteed to 
be real when a-\- 5 and 7 have the same phases, i.e., 

r = ie'^ 7 = ce''^. (9) 
Czech. J. Phys. 53 (2003) A 3 
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The real spectrum point interaction ((SJ is parameterized by 6 real parame- 
ters. The four-parameter family of self-adjoint (non-separated) boundary condi- 
tions is contained in this 6-parameter family. The family of PT-symmetric 
(non-separated) boundary conditions leading to operators with real spectrum is 
also included in the family of boundary conditions JSj) or @ . 



4 Integrable many-body systems 

The self-adjoint boundary conditions ^ and PT-symmetric boundary con- 
ditions © and Q), and the real spectrum boundary conditions Q with the parame- 
ters satisfying or Q also describe two spinless particles moving in one dimension 
with contact interaction when they meet (i.e. the relative coordinate x = 0). When 
the particles have spin s but without any spin coupling among the particles, ip rep- 
resents any one of the components of the wave function. In the following we study 
the integrability of one dimensional systems of A^-identical particles with general 
contact interactions described by the non-separated boundary conditions that are 
imposed on the relative coordinates of the particles. We first consider the case of 
two particles (N — 2) with coordinates xi, xi and momenta fci, ^2 respectively. 
Each particle has n-'spin' states designated by s\ and S2, \ < Si<n. For x\ ^ X2, 
these two particles are free. The wave functions ip are symmetric (resp. antisymmet- 
ric) with respect to the interchange (xi, si) {x2, S2) for bosons (resp. fermions). 
In the region xi < X2, from the Bethe ansatz the wave function is of the form, 

ip = Q,^2e*(''i="i+'=^^^) + a2ie*"='"''+'''"''\ (10) 

where ai2 and a2i are x 1 column matrices. In the region xi > X2, 

f = (P^^ai2)e*('^'l^^+'^^^l) + (Fl2Q,2^)g»(fe2a=2+fci:.l)^ (11) 

where according to the symmetry or antisymmetry conditions, P^^ — p^^ for bosons 
and P^^ = —p^^ for fermions, p^^ being the operator on the n? y. 1 column that 
interchanges si ^ S2. 

Set ki2 = {ki — fc2)/2. In the center of mass coordinate X = {xi + a;2)/2 and 
the relative coordinate a; = a;2 — we get, by substituting and 1)111) into the 
boundary conditions (|7|) at a; = 0. 

r ai2 -f a2i = a-P^^(ai2 + a2i) + i/3fci2^'^^(ai2 - Q;2i), ,.ry. 
\ iki2{a2i - 0-12) ^ "fP^'^{ai2 + a2i) + i5ki2P^'^{cti2 - ct2i)- 

Eliminating the term P^^ai2 from (|12ll we obtain the relation 

a2i = Yi^ai2 , (13) 



where 



12 _ 2iki2iaS - /^t)^ + ^fci2(« -S) + (fci2)^/3 + 7 

*fci2(a + 5)-H(fci2)2/3-7 ■ ^ ' 
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For > 3 and xi < X2 < ■■■ < xn , the wave function is given by 

— ^^piikixi + k2X2 + ... + kMXN) \ „i(fe2a;i + fei2;2 + ... + fejva;jv) 

+ {N\ - 2) other terms. ^ ' 

The columns a have x 1 dimensions. The wave functions in the other regions 
are determined from (|15|l by the requirement of symmetry (for bosons) or antisym- 
metry (for fermions). Along any plane Xi = Xi+i, i € 1,2, — 1, from similar 
considerations we have 

ahh...iih+i...iM = ^"^t/i"'i'2--i.+ii.--'jvi (16) 

where 

y^^+l _ 2»fcM.+l i^S - (3j)P" + ' + {a-d) + (fcM.+l ?P + 7 , . 

'•+^'* ikuu+Aa + S) + {kuu+,rP-l ■ ^ ' 

Here ki^i^^^ = {ki. — ki.^^)/2 play the role of spectral parameters. P"+i = 

for bosons and P"+i = — for fermions, with the operator on the x 1 

column that interchanges Si Si+i. 

For consistency Y must satisfy the Yang-Baxter equation with spectral param- 
eter IHHH, i.e., 

T^m,m+lTAm+l,m+2TAm,m+l ^ T^m+l,m+2TAm,m+lTAm+l,'^?i+2 

ij kj ki ki kj ij ' 

or 

ij kj ki ki kj ij \ " / 

if m, r, 5 are all unequal, and 

'^3 3'^ ' ijkt kl IJ V/ 

if m, r, s, q are all unequal. These Yang-Baxter relations are satisfied when 

aS-Pj^l, /? = 0, a = 6, (20) 

i.e., (3 = 0,a = S = ±l. 

Therefore for self-adjoint contact interactions Ql, the iV-body system is inte- 
grable when 6 = 0, a — d = ±1, b — 0, c arbitrary. The case a — d — l,9^b = 
corresponds to the usual (5-function interactions, which has been investigated in 
|S, J . The case a = d = —1, 9 = b = 0, is related to a kind of anti-(5 interactions 

m 

For iV-body systems with PT-symmetric contact interactions (O , the integrable 
condition (|20|1 implies that — (p ~ b — 0, which is just the usual self-adjoint S- 
interaction. 

From @ and H2U|I . we also have that an many-body system with contact inter- 
action and real spectra is integrable only when a — 6 — ±1, j — c for some c € JR, 
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Fig. 1. Relations among point interactions according to integrability. 

which is again the self-adjoint J-type interactions, see Fig.l for the relations among 
point interactions according to integrability. 

We have presented a complete picture of self-adjoint, PT-symmetric and real 
spectrum point interactions, and their corresponding integrability. What we con- 
cerned here are just the case of particles with only pure contact interactions, and 
the possible contact coupling of the spins of two particles are not taken into ac- 
count JS]. A further study along this direction would possibly give rise to more 
interesting integrable quantum many-body systems with various symmetries and 
spectrum properties. 
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